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We consider a double quantum dot system with two embedded and non-aligned spin impurities 
to manipulate the magnitude and polarization of the electron spin density. The device is attached 
to semi-infinite one-dimensional leads which are treated exactly. We provide a real-time description 
of the electron spin dynamics when a sequence of ultrafast voltage pulses acts on the device. The 
numerical simulations are carried out using a spin generalized and modified version of a recently 
proposed algorithm for the time propagation of open systems [Phys. Rev. B 72, 035308 (2005)]. 
Time-dependent spin accumulations and spin currents are calculated during the entire operating 
regime which includes spin injection and read-out processes. The full knowledge of the electron 
dynamics allows us to engineer the transient responses and improve the device performance. An 
approximate rate equation for the electron spin is also derived and used to discuss the numerical 
results. 

PACS numbers: 73.63.-b,72.25.-b,85.75.Mm 



I. INTRODUCTION 

The ability of controlling magnitude and orientation of 
electron spin densities in integrated molecules and quan- 
tum dots is of utmost importance to bring quantum com- 
putation closer to real lifeP^The microscopic description 
of nanoscale spin devices like, e.g., the two- quantum-bit 
gate envisaged by Loss and Di Vincenzop' constitutes a 
challenging problem in the theory of open systems far 
from a steady state. Research activities in the emerging 
field of spin-dependent transport^ 4 have mainly focussed 
on steady-state properties. Only very recently the tran- 
sient dynamics of spin polarized curre nts thro ugh quan- 
tum dots has attracted some attention] 5 * 6 * 7 * 8 ! partly due 
to experimental advances in manip ulating electronic den- 
sities with ultrafast voltage pulses.^^ 12 ' 13 ' 14 ' 15 ' This 
paper goes in the same direction and wants to be a further 
step toward the bridging of spin dependent transport and 
fundamental quantum computation. We perform time- 
dependent simulations of the charge and spin dynamics of 
a nanoscale device in contact with one-dimensional leads. 
The semi-infinite leads are treated exactly. The results 
are analyzed within the framework of non-equilibrium 
Green's functions. 

We consider a double quantum dot device to manip- 
ulate the spin orientation of spin-polarized electrons. 
Both quantum dots contain a static spin impurity with 
which the electron spin is coupled. The exchange cou- 
pling constant is much larger than experimentally ac- 
cessible Larmor frequencies, a feature that renders the 
spin impurity a potentially ultrafast mean to rotate the 
electron spinpiiLLZl Model systems of quantum transport 
through magnetic quantum dots have been previously 
used to study the conductance oscillations of a local 



nuclear spin in a magnetic field, 18 the gauge- invariant 
nature of the charge and spin conductances,^ the spin- 
interference and Aharonov-Bohm oscillations ina quan- 
tum ring with embedded magnetic impurities j 2Q l 21 l 22 -l and 
the effects of the entanglement of two spin impurities on 
the conductance! 23124 1 

In this paper we focus on the short time response of the 
system when subject to a sequence of voltage pulses, as 
illustrated in Fig. [I] The injection of spin polarized elec- 
trons from the left lead to the first quantum dot (QD1) 
is followed by a rotation of the electron spin in QD1. 
Afterwards the electron spin is transferred from QD1 to 
the second quantum dot (QD2) and its polarization is 
maintained parallel/antiparallel to the spin impurity of 
QD2. Eventually, the electron spin in QD2 is read out 
by calculating the spin current at the interface with the 
right lead. We provide a time-dependent description of 
some crucial processes in the theory of spin transport, 
namely the injection of spins from a lead to a quantum 
dot and the spin dynamics of a double quantum dot sys- 
tem weakly coupled to leads. The results of our analy- 
sis include 1) an overshooting of the spin accumulation 
during the spin injection phase, 2) a considerable delay 
in the spin relaxation for different exchange couplings in 
QD1 and QD2, and 3) oscillations in the transient regime 
whose frequency depends on the bandwidth of the leads 
and, therefore, are absent in the commonly used wide 
band limit approximation. 

The paper is organized as follows. In Section [II] we 
present our model system and introduce the basic nota- 
tion. A set of approximate equations to describe different 



operating processes are derived in Section [nil We obtain 
a rate equation for the electron spin of a quantum dot in 
contact with an electron reservoir and identify the mech- 
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FIG. 1: Schematic of the double quantum dot device coupled 
to leads. At t < To = the system is in equilibrium. At To a 
spin bias Ul,i is switched on and simultaneously the barrier 
between QD1 and the left lead is lowered. The injection of 
spin up electrons ends at T\ when the bias is turned off and 
the barrier is raised up. Now the spin in QD1 rotates till T2 
when the barrier between the dots is lowered and the electron 
spin is transfer from QD1 to QD2. At T3 the interdot barrier 
is raised up again while the barrier between QD2 and the 
right lead is lowered. Tuning the electrochemical potential in 
the right lead Ur to be in between the two levels of QD2, 
we measure a large (small) spin current if the electron spin is 
parallel (antiparallel) to the spin impurity. 



anisms leading to a deterioration of the spin polarization 
and to a damping of the spin magnitude. This analysis 
will then be used to optimize the spin injection from one 
of the leads to one of the quantum dots. We also investi- 
gate the spin transfer between the two quantum dots for 
different initial orientations of the electron spin. In See- 



the algorithm of Ref. [25j see Appendix [Bj to perform nu- 
merical simulations of the microscopic electron dynamics 
of the double quantum dot system. The results are then 
interpreted and discussed using the framework developed 
in Section [TTT1 In Section [V| we summarize the main find- 
ings and discuss future directions. 



II. MODEL SYSTEM 

We consider a basic two spin-impurities model consist- 
ing of two one-dimensional leads coupled to two single- 
level (per spin) quantum dots. The first quantum dot is 
connected to the left lead (L), the second quantum dot is 
connected to the right lead (R) and a hopping term ac- 
counts for tunneling of electrons between QD1 and QD2, 
see Fig. IT] The Hamiltonian describing the left (L) and 
right (R) leads is 



V C ]-\-la,a C ja,a 



*=U 3=0 



(1) 



with a = L,R. In Eq. ([I]) the quantity V is the hopping 
integral between nearest neighbors orbitals and e a ^ a (t) 
is the time-dependent on-site energy of lead a which, in 
general, can depend on spin. For e a ^ = the energy 
window of both L and R continua is (— 2|V|,2|V|) and 
the half-filled system correspond to a chemical potential 
/ihf = 0- The Hamiltonian of the double quantum dot 
system reads 



QD 



JiSi ■ ^ dt,a^cra'd it<7 ' + Vi(t) ^ d\^d i)0 



+ E [^QD(t)<^ + ^QD(*)4* d l."] 



(2) 



with Si = (sin Oi cos fa, sin Oi sin fa, cos Oi) the spin of the 
impurity i = 1, 2, Ji > the corresponding antiferromag- 
netic exchange coupling, and a = (a x , a y , a z ) the Pauli 
matrices. The tunnel barrier between the dots can be 
tuned by varying an external gate voltage. 3 This is mod- 
elled by using time-dependent gate voltages v\(t), V2(t) 
and interdot hopping integral Vqd(£)- The first term in 
Eq. ([2| is the Hamiltonian of the two isolated QDs and 
can conveniently be rewritten in matrix form as 



E«t>4i) 



Jie 



Ji cos Oi 
sin<9 ? - 



Jie*** sin Oi 
Vi — Ji cos Oi 



H,1 



tion IV we use a spin-generalized and modified version of 



(3) 

From Eq. (|3]) we see that the isolated QD has two levels 
at energy £^± = Vi =b Jj. If Ji > \vi\ and e aj(7 = one 
level is above /ihf while the other is below. 
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The double QD system is connected to the left and 
right leads via the tunneling Hamiltonian 

H T = ^[v L (t)4, <r c L,* + V£(t)4 Lta d 1 , <r '_ 

a 

+ J2 [V R (t)4 t<r C0R,a + vs(t)4 fli< A„] • (4) 



As for the interdot coupling Vqd we allow the hopping 
integrals and Vn(t) to be time-dependent. 

Below we discuss a sequence of operations to manipu- 
late the orientation of the electron spin in QD2 for a fixed 
orientation of the spin of the injected electrons. With- 
out loss of generality we assume that for negative times, 
t < 0, the whole system is in equilibrium at chemical 
potential \i and inverse temperature /3 [Fig(]Ji)]P^ The 
two quantum dots are initially very weakly coupled to 
the leads, V a « V, and between them, Vqd << V. 
Furthermore, the two energy levels of both QD1 and 
QD2 are much larger than the chemical potential, i.e., 
Si^± » /i, and the population on the dots is practi- 
cally zero. Starting from this configuration we apply a 
sequence of four perturbations to 1) inject spin up elec- 
trons on QD1 [Figflj})] 2) rotate the electron spin in QD1 
[Fig 1 z)] 3) transfer the electron spin from QD1 to QD2 
[Fig 1 i)] and 4) read out the polarization of the electron 
spin in QD2 [Fig[]J)]. Due to the wide range of possi- 
ble time-dependent perturbations we restrict the analysis 
to piece-wise constant (in time) parameters and obtain 
a set of approximate equations to study the four differ- 
ent processes. This study will then help us in selecting 
the parameters for target-specific numerical calculations. 
Full simulations of the entire sequence will be shown in 
Section ED 



tion for the lesser Green function , on QD1 
. d 



i-G < (t-t) = [H Q v ll G < (t-t)}- 
(t; t)G A (i; t) + E£(t; t)G < (f; t) + h.c. 



dt 



(5) 



where we use boldface to indicate 2x2 matrices in spin 
space and the symbol "[,]" denotes a commutator. In 
the above equation -H"qdi = v\ + J\S\ • a is the one- 
particle Hamiltonian of the isolated QD1 while Xlz, is the 
embedding self-energy of lead L. We have discarded the 
integral between and —i/3 along the imaginary time 
axis since V L (t < 0) = V (0) « V and hence E L ~ 
in equilibrium. 28 The superscripts R/A in 5]^ and G de- 
note retarded/ advanced components. The self-energy is 
diagonal in spin space since there is no spin-flip hopping 
between lead L and QD1. In terms of one-particle eigen- 
states ipk(j) and eigenenergies Sk of lead L one finds for 
M' > 



du _■ 



27T 



l^(Q)l 2 



(6) 
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27r^|^ fe (0)| 2 <5(w-e fc ). 



(7) 



At low temperatures and low biases only frequencies close 
to the Fermi energy e F are probed. Using the Wide Band 
Limit (WBL) approximation, i.e., 



^|V^(o)| 2 (^-^)~p F 



(8) 



III. THEORETICAL FRAMEWORK 

A. Spin injection and spin read-out 

At time t = we inject spin- up electrons into QD1 by 
suddenly switching on a spin bias, 27 SL,a(t) — 0(i)UL,a^ 
and reducing the height of the barrier between lead L 
and QD1, i.e., V L (t) = 6(-t)V^ + 6{t)V^ . The in- 
jection process terminates by switching off the spin bias 
and raising up the barrier to the equilibrium value V^ ^. 
There are two different mechanisms which contaminate 
the spin-up injection with x and y components. The first 
is the spin-precession around the spin impurity Si while 
the second is the spin-relaxation due to the increased 
electron hopping Vl- To tackle this problem we take ad- 
vantage of the fact that QD1 and QD2 are, during this 
process, weakly linked and we only consider the electron 
dynamics on QD1 in contact with the left reservoir, i.e., 
we approximate Vqd = 0. Using the non-equilibrium 
Green's function formalism one finds the following equa- 



with pY = p(sf) the local density of states at the inter- 
face, we can approximate Eqs. (plFFl as 



(9) 



^ aa ,{t-t') = --5 ai7 ,T5{t-t'), 



Z< aa ,(t;t') = i6 acr ,T J 



2tt 



e -i{u+U L ,„){t-t f ) f ^ 



Z7T 



-i£F,cr(t — O 



--i*5(t-if) 



(10) 



with T = 2tt\VW \ 2 p F and e F , a = e F + U Ly(7 . In Eq. (To} 
we have further approximated X)^ with its expression at 
zero temperature. 

Inserting these results into Eq. ([5| one obtains 



^G<(t;t) = [H QD1 ,G<(t;t)} - % -{T,G<(t;t)} 



2tt Jo 



* exp[-iS F (t - t)] A 



dt 



t - 1 



G A (t;t) + h.c. 



.(11) 
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where the symbol "{, }" denotes the anticommutator and 
the matrices [T] a(T > = 5 aa >T and [£f]<t<t' = (W^f,^ 
From Eq. (11) we can extract a rate equation for the 
electron spin 



52, el during the read-out phase is identical to Eq. ( 14) for 



Si, ei even though the parameters are different and, more 
importantly, different initial conditions must be imposed. 



S 1M (t) = --Tr[G<(t;t)a] 



(12) 



on QDl. In the WBF approximation the advanced 
Green's function reads 



G A (i;t) = iO(-At)exp[-i(H QD1 + -T)At] 



cos( Ji At) - i sin( J x At) Si • 3 , (13) 



with At = t — t. Substituting this result into Eq. (11), 



multiplying with a and tracing over the spin indices we 
find 

As,,, = .MS, as,,,)- ^s*,,., 

x [cos(e_At)sin(JiAt)6 r i 
- sin(6:_At){cos(JiAt)i - sin(JiAt)i A Si}], 

(14) 

with z the unit vector in the z direction, and e± = £F,|i 
£f,i- It is instructive to expand the right hand side in 
powers of t. To first order in t one finds a simplified rate 
equation for the electron spin 



dt 



Sl, e 



Jl(5i A Si, el) — ~^Sl,el 



7T V 



Te-z- Jir5i 



which is reliable for times t 1 > max[e_ 



(15) 

,Ji]/27r. From 



the above equation we can identify four different contri- 
butions. The term proportional to z is the spin-injection 
term and is responsible for an increase of the electron 
spin along the z direction. Such increase is quadratic in 
time and faster the larger the difference S- = Ul,\ — Ul,i 
is. The first and the last terms are responsible for a de- 
terioration of the spin direction due to spin precession 
(first term) and spin relaxation (last term). The latter 
drives the electron spin towards a configuration antipar- 
allel to the spin impurity S±. Finally, the second term is 
responsible for an overall damping of the spin magnitude. 
Going beyond the first order in t, see Eq. (14), one ob- 
serves the appearance of a new relaxation direction, that 
is z A Si. This latter result is completely general as it is 
only dictated by the symmetry of the system. 



We wish to emphasize that the rate equation ( 14 ) has 



been derived under the sole assumption that the quantum 
dot QDl is initially isolated and then contacted with lead 
L. This is the same situation occurring in the spin read- 
out phase when the barrier between the weakly coupled 
QD2 and lead R is lowered. Thus, the rate equation for 



B. Spin rotation and spin transfer 

After a time Ti the spin bias is switched off and the 
hopping Vl is again reduced to values much smaller than 
V. In this phase QDl is well isolated and the electron 
spin precesses around the spin impurity Si according to 



dt 



Si, el = Jl(Si A Si, el), t > Ti. 



(16) 



Let us now specialize to the situation illustrated in 
Fig. [l] with Si oriented along the positive x axis and 
§2 along the positive z axis. We recall that the Fermi 
energy is much smaller than the energy levels of the two 
isolated quantum dots and hence that the equilibrium 
electron density is vanishingly small. For Ji » e_, 
see Eq. (15), we expect an efficient injection of spin up 



electrons in QDl and for Ji »Ta major contamination 
along the y direction. This implies that the electron spin 
£i,ei(^i) nas a small x component at the end of the spin 
injection process. Since Si is parallel to the x axis Si, e i(t) 
rotates in the yz plane for t > Ti. We let the system 
evolve till a time T 2 > Ti and we approximate Si^iO^) = 
(0,5? >el (T 2 ),5f >el (T 2 )) on the yz plane and S 2 J(T 2 ) = 
(this latter approximation comes from the fact that Vqd 
and Vr are both much smaller than V for t < T 2 ). 

At t = T 2 we transfer the electron spin by lowering 
the barrier between QDl and QD2. This corresponds 
to an increase of the interdot hopping Vqd- Letting 
|$(T 2 )) be the evolved many-particle state of the en- 
tire system at t = T 2 , the density matrix p of the dou- 
ble quantum dot system has matrix elements [p]ia,ja' — 
($(T 2 )|dt ,d ijCr |$(T 2 )). It is convenient to introduce the 
notation 1 = (1,|), 2 = (1,|), 3 = (2, |) and 4 = (2, |) 
for the collective index (z, a). The density matrix is then 
represented by the following 4x4 matrix 



/ 



p = 2Si,el(T 2 ) 



with sin 20 



cos u 
| sin 20 







-f sin 20 \ 



sin 2 







0/ 



(17) 



S?el( T 2)/Sl,el(T 2 ), COS20 



Sf >el (T 2 )/Si >e i(T 2 ), and 5i, e i = yf Si M • S hel the 
spin magnitude. We are interested in how to choose the 
angle in order to maximize the electron spin of QD2 
along the z direction (parallel/antiparallel to the spin 
impurity 5 2 ). For simplicity we take the gate voltages 
vi = and v 2 = 0. Then the isolated double quantum 
dot system is described by the 4x4 Hamiltonian matrix 



H 



QD 



/ Ji(7 x Vqd1 2 
V VqdI-2 J2°z 



(18) 
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with I2 the 2x2 identity matrix. In terms of p and 
Hqv the z component Sf e i °f ^he electron spin on QD2 
is given by 

Sf,ei(i + T 2 ) = ^Tr [S| exp(iilQ D t)pexp(-iJI QD t)] , 

(19) 

with the spin operator of QD2 



2 2 

O2 cr z 



(20) 



and O2 the 2x2 null matrix. Substituting p from Eq. 
((171) we nnd 



S1AT2) 



= -sin2fl{[Si(t)] 1 , 2 -[S|(t)] 2)1 } 



2 

cos 2 0[£|(t)] 



1,1 



•sin 2 0[E|(t)] 2 , 2) 



(21) 

where we have defined the spin operator in the Heisen- 
berg representation 

E£(t) = exp(-iff QD *)£2 exp(zJJ QD *). (22) 
It is easy to prove that the function 0(t) = |(p2(^)]i,2 — 
[5^2 (^)] 2,1) is an °dd function of time while [E2(£)]i,i and 
[S2(t)]2,2 are even functions of time. In appendix |A| we 
further prove that 

E(t) = \E* 2 (t)h,i = -[S|(t)] 2 ,2, (23) 

which leads to the simple formula 

Slei(t + T 2 ) = S he i(T 2 )[0(t)sm26 + E(t)cos26] 

= 0(t)Sl el + E(t)Sl e] . (24) 

The function E(t) can be written as a linear combi- 
nation of the cosine functions cos^^t) while 0(t) as a 
linear combination of the sine functions smiu^t), where 
^iiv — Zfi — Zv is the difference between two eigenvalues of 
Hqe>. The eigenvalues £ M , \i = 1, . . . , 4 can be calculated 
analytically and read 



\ 



(25) 



with J\ = J\ ± J 2 2 . 

As an example, in Fig. |2j we plot the ratio 5| e i(t + 
T2) / Si^ e \(T2) as a function of time £ and initial polariza- 
tion for Ji = J2 =0.1 and Fqd = 0.2. We notice 
that for most polarizations S% G \(t + T2) remains smaller 
than 0.2. Only for some special value of the z compo- 
nent of the spin in QD2 reaches a value larger than 0.4. 
This means that the maximum efficiency in transferring 
an electron spin polarized in the yz plane from QDl to 
QD2 with final polarization along the positive z axis is 
about 80 - 90%. 

All above processes can be numerically simulated with- 
out resorting to any of the approximations employed in 
this Section. This allows for a more quantitative investi- 
gation of the device performance and will be the topic of 
the next Section. 




FIG. 2: Ratio S^ e i(t + T 2 ) / Si, e \(T 2 ) versus time t and initial 
polarization for J\ — J 2 — 0.1 and Vqd = 0.2. 



IV. NUMERICAL SIMULATIONS AND 
DISCUSSION 



In this Section we use a modified version of the algo- 
rithm proposed in Ref. [25] to propagate in time finite 
systems in contact with infinitely long leads, see Ap- 
pendix [Bj and investigate the microscopic dynamics of 
the spin-injection, the spin-accumulation as well as the 
spin rotation of conducting electrons scattering against 
the double QD device of Eq. ([2|. In the following 
analysis energies are measured in units of V, times in 
units of h/V, spins in units of H and currents in units 
of eV/h, with e the electron charge. The full Hamil- 
tonian is time independent for negative times and the 
system is in equilibrium at zero temperature and Fermi 
energy £p- We start by considering two identical QDs 
with exchange coupling J\ = J 2 = 0.1 and gate poten- 
tial v\ = V2 = weakly coupled to the left and right 
leads, V L = V R = = 0.01, and with interdot hop- 
ping Vqd = 0.01. Choosing V ~ 100 meV the exchange 
couplings Ji, J2 ~ 10 meV lie in the physical parameter 
range 17 and the corresponding time unit is h/V ~ 100 
fs which is appropriate to study ultrafast dynamics! 29 1 30 1 
The impurity spin Si of QDl is oriented along the posi- 
tive x axis while §2 is oriented along the positive z axis. 
The on-site energies of the leads e ajCr are initially all zero. 



A. Spin injection 

At time t = we switch on a spin bias £l,tW = 
0(£)^£,T m l ea d L for spin up electrons (Ul,i = 0) and 
increase the hopping Vl from V^ to V^\ 

In Fig. [3] we study the spin- inject ion process when 
the Fermi energy is Sp = —0.96 (which correspond to an 
initial electron occupation on QDl of the order of 10 -5 ) 
and the hopping between L and QDl at positive times is 
= 0.5. We calculate the time-dependent expectation 
value of the spin of the conducting electrons Si je i on QDl 
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FIG. 3: The three components of the electron spin in QD1 
versus time for a sudden change of the hopping Vl from 
y(°) = 0.01 to — 0.5 and a simultaneous sudden switch- 
ing of the spin-bias Ul,i =0.7, 0.9, 1.1, 1.3. The equilibrium 
parameters are J\ — J2 = 0.1, Vl — Vr — Vqd = V^ and 
the Fermi energy is ef — 0.96. 



for different biases Ul^- Since Ul^ ~ 1 > J\ the rate 
equation (15) is reliable for times t < 2tt/Ul^ ~ 27r. 



In this time window we observe that the z component 
Sf el increases quadratically in time and that the rate 
is larger the larger is the spin bias Ul,^ in agreement 
with Eq. (15) (we recall that in this case S- = Ul^)- 
The y component Sf el has a trend similar to 5f el but 
the transient is even smoother. This can be explained 
by observing that as the spin up electrons enter QD1 
they undergo a spin rotation due to the spin impurity 
oriented along the positive x axis. Taking into account 
that for small t we have S* el ~ t 2 , from Eq. (15) we 
see that Sf el ~ Jit 3 . As the z component also the x 
component Sf el grows quadratically in time. From Eq. 
(15) one finds Sf el (t)/Sf el (t) ~ S-/J\ meaning that to 
minimize the contamination of spin up electrons with an 
x component it must be S-j J\ > 1. 

We wish to observe that at intermediate biases Ula 



the numerical results agree with the rate equation (14) 
only qualitatively. The comparison between the time evo- 
lution of the electron spin in QD1 for Ul^ = 0.7, 0.9 
as obtained with one-dimensional leads and with leads 
treated in the WBL approximation is shown in Fig. [4j 

In Fig. [5] we fix the bias for spin-up electrons to be 
Ul,i = 1-3 and analyze the spin-dynamics on QD1 for 
different values of V^\ We first observe that the tran- 
sient time decreases by increasing V^ and hence T. This 
is easily understood by noticing that the second term in 
Eq. (14) yields an exponential damping of the spin os- 
cillations. The spin oscillations can be observed in the 
y and z components and are due to the spin precession 
around the spin impurity S±. The period of the oscil- 
lation is T = tt/Ji and is independent of V^\ as it 
should. It is also interesting to observe that for small 
times S* e i overshoots its steady value and hence more 



FIG. 4: Comparison between the results obtained with one- 
dimensional (ID) leads and WBL leads for the electron spin 
in QD1 and Ul^ — 0.7, 0.9. The other parameters are the 
same as in Fig. [3] 
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FIG. 5: The three components of the electron spin in QD1 
versus time for a sudden switching of the spin-bias Ul,i — 1.3 
and a simultaneous sudden change of the hopping Vl from 
= 0.01 to = 0.1, 0.2, 0.3, 0.4. The other parame- 
ters are the same as in Fig. [3] 



efficient spin injections may be achieved by properly en- 
gineering the transient response. In our case, for an effi- 
cient spin up injection only the ratios r xy — \Sf e \/ S\ el | 
and r xz = \Sf e \/S* el | must be small at the end of the 
process since the y component can be reduced to zero 
in the second phase when Vl « V and 5i je i can pre- 
cess around the spin impurity. From Fig. [5] we find that 



for V par = 0.2 and at t 



10 the ratios r xy ~ 0.1 and 



r X z ~ 0.12 while at the steady state r xy ~ 0.31 and 
r xz - 1.43. 



B. Spin rotation 

The injection process terminates after some time T\ by 
switching the spin bias off and raising back the barrier 
between QD1 and the left lead, i.e., V L (t) = 0(-t)V w + 



7 



1 1 1 


1 1 1 
a)" 


1 


J - x - 

- y 

z 

1 i T 




time 



FIG. 6: Electron spin in QD1 and electron density at the left 
interface during the injection phase (t < Ti) and the rotation 
phase (t > Ti) when Ti = 10. The equilibrium parameters are 
the same as in Fig. [5] i.e., Ji = J2 = 0.1, Vl = Vr = Vqd — 
y (0) and the Fermi energy is e F = 0.96. For < t < Ti 
the system is perturbed by a spin bias Ul^ — 1.3 and a 
larger hopping Vl = V^ 1 ) = 0.2. At t = Ti the spin bias 
is switched off and the hopping Vl is suddenly changed to 
y( 2 ) = 0.01. Panel a) The three components of Si, e \- Panel 
b) The trajectory of the projection of Si >e i onto the yz plane. 
Panel c) Spin up and down density on the first site of the left 
lead. The inset is a magnification of both densities ul^ and 
riLi for times t > 25. 



0(t)0(Ti -t) V^ +0{t-T{) V (2) . During the second phase 
QD1 is weakly coupled to the environment and the elec- 
tron spin precesses around S\ . Let us focus on the situa- 
tion discussed above with = 0.2 and let T\ = 10 be 
the duration of the first phase. In Fig. [6] we study the 
electron spin on QD1 [panels a) and b)] and the densities 
riL,a = (cq L a coL,a) on the first site of the left electrode 
[panel c)] for = 0.01. The contaminating component 
Sf el ceases to decrease at t = Ti while the y and z com- 
ponents are well described by damped cosine functions 
with a phase lag of 7r/(4Ji) [panel a)]. Due to the weak 

contact the magnitude 5i ?e i = \J *5i, e i * *5i,ei of the 
electron spin changes on a time scale much longer than 
the spin-exchange time-scale ~ 1/Ji- This is shown in 
panel b) where the trajectory of Si >e i is projected onto 
the yz plane. For times t < T\ the trajectory has a large 
radial component while for t > T\ the spin moves along 
a spiral trajectory. It is also interesting to look at the 
densities on the nearest neighbor site of QD1 [panel c)]. 
During the first phase (t < T\) a majority of spin up 
electrons are transferred from lead L to QD1 and, as a 
consequence, n L ^ decreases. On the contrary the density 
ul,i increases due to the following two-step mechanism. 
As the spin up electrons hop from lead L to QD1 they 
undergo a spin rotation and acquire a down component. 
These electrons have about zero energy and can easily 
hop to the left lead where the spin-down band is filled 
up to €f + Ul,i = £f = —0.96. At the end of the injec- 



tion process the densities change abruptly and approach 
their initial value since 

y{2) = y(0) The inget of panel 

c) is a magnification of the curves riL^(t) and n^^if) for 
25 < t < 160. It is clearly visible a quantum beating 
in both densities due to the alignment of the spin impu- 
rity along the x axis. In both cases two oscillations with 
frequency \ep ± J\ \ = 0.96 ± 0.1 are superimposed to an 
envelope oscillation of frequency 2Ji = 0.2. 




FIG. 7: The three components of S\^\ [panel a)] and the 
trajectory of the projection of Si je i onto the yz plane [panel 
b)] for the same system as in Fig. [6] except that the hopping 
parameter — 0.06 is six times larger. Panels c)-d) like 
panels a)-b) but the equilibrium parameter J2 = 0.02 while 
the hopping parameter — 0.01. 



The spin rotation phase is further investigated in Fig. 
[7] where we consider the same system as in Fig. [6] except 
for the value of the hopping parameter = 0.06 which 
is six times larger [panels a) and b)] or the exchange cou- 
pling J r 2 = 0.02 which is five times smaller [panels c) and 
d)]. In the first case the x component remains an order 
of magnitude smaller than 5i ?e i [see panel a)] and even- 
tually approaches a steady value slightly larger than the 
initial one [not shown]. As in Fig. [6] the electron spin is 
damped in all three directions but it decays faster. The 
projection of 5i ?e i onto the yz plane [panel b)] yields a 
spiral trajectory which finishes very close to the origin 
after a time t ~ 160. On the other hand, for a smaller 
coupling J2 = Ji/5 we do not appreciate any damping 
within the time propagation window t < 160. The y 
and z components of Si je i are well described by two un- 
damped cosine functions with a phase lag 7r/(4Ji) and an 
amplitude which is about ten times larger than \Sf el |, see 

panel c). In panel d) we show the projection of *5i ?e i onto 
the yz plane. The reduced damping is a desirable feature 
and has to be attributed to the mismatch of the energy 
levels in the two quantum dots: ± J\ in QD1 and ± J 2 in 
QD2. 
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FIG. 8: The electron spin S 2 , e \ in QD2 before (t < T 2 ) and af- 
ter (t > T 2 ) the spin transfer phase when the interdot hopping 

0.5 [panels 



is set to Vqq =0.2 [panels a) to c)] and V c 



(i) 

QD 



d) to f)]. The transfer phase starts at T2 ~ 36.5 [panels a) 
and d)], T 2 ~ 40.4 [panels b) and e)], and T 2 ~ 44.3 [panels 
c) and f)]. All parameters before time T2 are the same as in 
Fig. [6] except for the exchange coupling J 2 = 0.05. 



C. Spin transfer 

The rotation of the electron spin in QD1 (t > Xi) ter- 
minates at t = T 2 > T\ when the barrier between QD1 
and QD2 is lowered and, as a consequence, the interdot 



hopping increases, i.e., Vq D = V^0(T 2 - t) + V^6(t - 
T 2 ), where V$ = °- 01 - This is the 

spin transfer phase. 
In Fig. [8] we plot the three components of the electron 
spin in QD2 versus time for Vqq =0.2 [panels a) to c)] 
and V$ = 0-5 [panels d) to f)]. For times t < T 2 the 



system undergoes the same perturbations as in Figs. [6][7| 
Here we have considered an exchange coupling in QD2 of 
J 2 = 0.05 and a hopping = 0.01. The frequency 
of the oscillations is larger the larger the interdot cou- 
pling is, in agreement with Eq. (25). The efficiency of 



the transfer has been investigated for different times T 2 
at which the electron spin in QD1 is polarized along z 
[T 2 - 36.5, panels a) and d)], ^(z - y) [T 2 - 40.4, pan- 
els b) and e)], and — y [T 2 ~ 44.3, panels c) and f)]. For 
our choice of parameters the efficiency is higher if the 
spin in QD1 is polarized along z. 

We also observe that for all three components the max- 
ima of the electron spin in QD1 correspond to the minima 
of the electron spin in QD2, see Fig. N where we plot Si iG \ 
and S2,ei for T 2 = 36.5 and Vq^ =0.2 [panels a) to c)] 
and V$ = 0.5 [panels d) to f)]. From Fig. ^ panel d) 
and Fig. [9] panels d) to f) one observes that when T 2 
corresponds to the time at which *Si >e i (^2 ) is polarized 
along z, the maxima of 5| el are close to the zeros of 
S% e i and S% el , in agreement with the analysis of Section 



IIIB 



We define the ratio r_\_ = S< 



!,el/ \J (^2,el) 2 + (^2,el)^ 



In the propagation window the maxima of 5| el occur at 
t = 39.32 when r± ~ 0.28 and 5| el = 0.153,' t = 95.68 



when r± ~ 0.26 and S, 



FIG. 9: Electron spin in QD1 and QD2 for T 2 = 36.5 and 
V$ = O- 2 [panels a) to c)] and V$ = °- 5 [panels d) to f)]. 
The other parameters are the same as in Fig. [8] 



r± 



0.24 and 5; 



2,el 



0.151. Taking into account that 
Sl eX (T2) = 0.163 the efficiency of the spin transfer can 
be up to 90%. 



D. Spin read out 

At a time t = T3 when S^ei has a maximum or a 
minimum, the interdot hopping is lowered, i.e., Vqd(^) = 



V<$0(T 2 -t) + V$0{t - T 2 )8(T 3 -t) + V$0(t - T 3 ) with 

(2) 

Vqy) << V, and the spin transfer phase ends. 

In Fig. 10 we consider the same system parameters 
and perturbations of Fig. [9] (with Vqq = 0.5) and fix 
the time T3 = 39.32 when S% e \ has a maximum. At 

(2) 

t = T3 the interdot hopping is lowered to = 0.001 
and QD2 becomes an almost isolated system. At this 
stage the density of spin up and down electrons in QD2 
is practically constant as one can see from the insets 
of Fig. 10 in panels a) and b). Shortly after T3 the 
read out phase starts. At t = 60 we lower the bar- 
rier between QD2 and lead R and simultaneously switch 
on a bias Ur^ = Ur^ = Ur = 0.96 in the right 
lead. The electrochemical potential in lead R becomes 
I^r = £f + Ur = and lies in between the two energy 
levels £2,± = =t^2 = ±0.05 of the isolated QD2, with the 
highest level £ 2>+ for spin up electrons and the lowest level 
E2- for spin down electrons P- 1 Spin up electrons in QD2 
have, therefore, energy larger than \ir and tunnel to the 
lead R. As a consequence the spin up density decreases, 
as one can see in Fig. 



(2)/ 



10 panel a). On the contrary, 



the lowest level 82,- has energy below fiR and a vanish- 
ingly small occupation. Spin down electrons tunnel from 
lead R to QD2 and the density of spin down electrons 
increases, see Fig. 10 panel b). This charge transfer gen- 
erates a right-going spin- up current Ir^ and a left-going 
spin-down current Ir^, see Fig. 10 panel c), which re- 



sults in a large spin-current. The spin dynamics in the 
xy plane is displayed in Fig. 10 panel d) where, besides 



2, el 



0.152, and t = 151.96 when the monotonically decreasing z component, we plot the 



9 




100 200 300 400 500 

time 



100 200 300 400 500 600 

time 



FIG. 10: Spin up density [panel a)] and spin down density 
[panel b)] in QDl and QD2 (the inset shows a magnification 
of the time window < t < 60). The spin polarized current 
at the right interface, /r,<x(£), is displayed in panel c) in units 
of 10~ 3 , while panel d) shows S 2 ,e\(t). For t < T 3 = 39.32 
the system parameters and the perturbations are the same as 



Fig. [ijwi 



with V c 



(i) 



QD 



is lowered to the value Vq^ 



0.5. For t > T3 the interdot hopping 

0.001 and QD2 becomes a 
well isolated system. At t = 60 > T3 the hopping Vr(£) 
between QD2 and lead R is raised up to — 0.05 and 

simultaneously a bias Ur^ — Ur,i = 0.96 is switched on in 
the right lead. At this time the electrochemical potential in 
lead R is /jlr = sf + 0.96 = and lies in between the two spin 
levels of QD2. 



x and y components of £2, el- Due to the symmetry of 
the problem 5f el and 5| el oscillate around zero with an 
exponentially decreasing amplitude. 
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FIG. 11: Spin up density [panel a)] spin down density [panel 
b)] on QDl and QD2, spin polarized current at the right in- 
terface, Ir,o-j in units of 10 -3 [panel c)], and S2, e \ [panel d)]. 
The insets of panels a)-b) show a magnification of the density 
in the time window < t < 60. Same system parameters and 
perturbations of Fig. [I5|but T 2 = 52.05 and T 3 = 54.84. 

The situation corresponding to the antiparallel config- 
uration in QD2 is analyzed in Fig. 11 The difference 



with the previous case is that we let the electron spin 
in QDl rotate till is polarized along the negative z axis. 
The first time Sf el is minimum occurs at T2 = 52.05, 
see insets in panels a) and b). The spin transfer phase 
ends at T 3 = 54.84 with an efficiency of about 90%. This 
can be seen in the inset of panel b) where the spin down 
density of QD2 swaps with that of QDl in the time win- 
dow (T2,Ts). At t = T3 the system undergoes the same 
perturbations considered in Fig. IlO| Being the spin up 



level of QD2 scarcely populated the change in the spin up 
density [panel a)] and spin up current at the right inter- 
face [panel c)] is very small as compared to the parallel 
configuration. A small change is observed for the spin 
down quantities as well due to a population of about 0.3 
in the spin down level of QD2. Contrary to the paral- 
lel configuration set up, the z component S| e \ is negative 
when the read out phase starts and does not change sign, 
see panel d). 




FIG. 12: Discrete Fourier transform of the spin current 
/spin = Ir,i — Ir,i at the right interface for the parallel and 
antiparallel configurations of Fig. [lO]and Fig. |TT| respectively. 
The inset shows / sp m in time domain, with t in the range (65, 
640), which has been used to perform the discrete Fourier 
transform. Both / sp i n (a;) and I sp in(t) are in units of 10 -3 . 

The spin current I sp in(t) = lR,i(i) — lR,i(t) during the 
read-out phase (t > 60) is displayed in the inset of Fig. 12 
for the parallel and antiparallel configurations analyzed 
in Figs. IOJTT One observes an exponential decay with 
superimposed oscillations of frequency \e-p + Ur ± J2 1 = 
0.05, as expected. However, a closer inspection reveals a 
richer structure. In Fig. 12 we show the discrete Fourier 



transform of I sp i n (t) with t in the range (65,640). Be- 
sides the peak at uo = 0.05 there exist an extra peak at 
frequency uo = |2 — ev\ ~ 2.96 and an asymmetric peak 
structure at frequency uo = \sv + 2| ~ 1.04. The extra 
transient frequencies are due to the finite bandwidth of 
the leads since the energies +2 and -2 (in units of V) 
correspond the top and the bottom of the right band. 

In conclusion, we have shown how to propagate in time 
a spinful open quantum system subject to arbitrary time- 
and spin-dependent perturbations. The semi-infinite na- 
ture of the leads has been exactly accounted for. Full 
simulations of the microscopic charge and spin transient 
dynamics of a double quantum dot in its operating regime 
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FIG. 13: Spin current at the left and right interfaces for the 
parallel and antiparallel configurations of Fig. ^] and Fig. 
[TT] respectively. For the parallel configuration T = 0, T\ — 
10.0, T 2 = 36.5, and T 3 = 39.32, while for the antiparallel 
configuration T = 0, Ti = 10.0, T 2 = 52.05, and T 3 = 54.84. 



have been presented. Figure 13 summarizes how the de- 
vice works by displaying the spin currents at the left 
and right interfaces during the entire sequence of volt- 
age pulses. Different processing of the injected spin up 
current results in different spin currents at the right in- 
terface. 



V. SUMMARY AND OUTLOOK 

In the last few years we have witnessed an increas- 
ing interest on transient responses in quantum trans- 
port mainly due to their potential relevancy in molec- 
ular electronics, a field where molecular devices will pos- 
sibly operate under non- steady- state conditions. The 
main difficulty in the study of the short-time response 
of open quantum systems stems from the macroscopic 
size of the leads. Several approaches have been pro- 
posed to tackle this problem. Treating the leads in the 
WBL approximation allows for obtaining a simple i nte- 
gral equation for currents, densities, etc.] 32 * 33 * 34 * 35 * but 
lacks retardation effects. One-dimensional leads have 
been ap proxi mately treated within a Wigner-function 
approach 36 ^ 37 ! or by i ncluding only a finite number of 
lead unit-cells! 38 * 39 * 40 * Only recently it became possible 
to deal with the semi-infinite nature of the leads using a 
scheme based on wave-functions propagation 2 ^ or, alter- 
natively, other algorithms based on so lving the D yson- 
Keldysh equations in the time-domain. 41 42 43 44 l 45 l Few 
attempts to include electro n-cor relatio n 1 46 1 47 1 as well as 
electron-nuclear interaction d 48 * 49 * in the transient regime 
have also been made. 

In this work we have used a modified version of the 
propagation algorithm of Ref. 25, see Appendix [Bj and 
generalized it to include the spin degrees of freedom. We 
have proposed a double quantum dot system to manip- 
ulate the charge and the spin of the electrons. Numeri- 
cal simulations of the entire operating regime have been 



provided. These include some of the crucial steps in the 
theory of quantum computation, like, e.g., the injection 
of spins and their read out. 

The transient electron dynamics when a device is per- 
turbed by ultrafast voltage pulses is not only relevant 
to our microscopic understanding but an exploitable fea- 
ture to improve the device performance. This has been 



explicitly shown in Section IV the efficiency of the spin 
injection can be much higher during the transient than 
at the steady state. We also have found that for a given 
height of the barriers between lead L and QD1, QD1 and 
QD2, and QD2 and lead R, the damping of the spin mag- 
nitude during the rotation phase is much smaller for dif- 
ferent exchange couplings, i.e., J\ ^ J2, than for J\ = J2. 
This means that the spin relaxation can be substantially 
delayed using different quantum dots. 

Using the non-equilibrium Green's function formal- 
ism in the WBL approximation we have obtained a rate 
equation for both the spin-injection and spin read-out 
processes. For short times the rate equation becomes 
remarkably transparent and permits us to identify the 
mechanisms leading to a relaxation of the spin magni- 
tude and to a deterioration of the spin polarization. Go- 
ing beyond the WBL approximation results in a richer 
structure of the transient responses, as transitions be- 
tween the Fermi energy and the bottom/top of the band 
occur as well. 



As shown in Section |IV[ the possibility of simulating 
operational sequences like, e.g., that of Fig. [I] allows 
for a real-time study of fundamental processes not ac- 
cessible otherwise. Much more work is, however, needed 
before a systematic comparison with experimental data 
can be made. Accounting for intradot and possibly inter- 
dot electron-electron interactions is of crucial importance 
for describing, e.g., the Coulomb blockade or the Kondo 
regimes. The complications here stem from the neces- 
sity of including electron correlations in a time-dependent 
conserving manner, a progress which can be made either 
within the framework of many body theory 50 or within 
one-particle framewo rks li ke, e.g., time-dependent den- 
sity functional theory.^^ Developments in this direction 
have been made in steady-state situations by treating the 
correlation at the GW level .HnHHHHU 

Another fundamental issue to be pursued is the exten- 
sion to three-dimensional leads. This would allow for a 
proper treatment of the long-range Coulomb potential as 
well as for a realistic description of the atomistic struc- 
ture of the tunneling barriers. 

Finally, the recent experimental advances in attaching 
quantum dots to superconducting leads^ prompt for a 
generalization of the propagation algorithm to leads de- 
scribed by, e.g., BCS-like models. Such development will 
give us access to a completely new phenomenology due 
to the competition between the pairing interaction and 
the spin- flip interactions, a topic not yet explored in the 
transient regime. 
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The result in Eq. ( 23 ) is a consequence of the relative 



orientation of the spin impurity Si with respect to S2. 
By definition the quantity [^Wliji * s ^ ne ma trix 
element of the product of three matrices 



[E£(t)]i,i = 



J\Vx Vqd12 
Vqd12 JlVz 



o 2 2 

02 CFz 



Vqd12 Ji^z 



(Al) 



1,1 



Consider the unitary operator U = U g U z U x which consists of a rotation of both spin impurities around the x axis 

exp[— in a x / 2) O2 
2 

2 

2 exp[-i7ra z /2} 



by an angle 7r, U a 
jj _ 1 exp[-z7rcr 2 /2] 



r . /oi J 5 followed by a rotation around the z axis by an angle 7r, 

exp[ 17T(7 X / Z\ J 

, followed by a gauge transformation which changes the sign of the fermion 



/ i 2 Oo \ + 
operators on QD2, U g = y q 1 j ' ^ nsert ^ ons °^ ^he identity matrix U J U in Eq. (Al) gives 

P2(*)]l,l = 



-it I JlCJx V ® d1<2 \ , , it ( Jl ° x T/ QD 1 2 j 

U ] Ue \Vqd12 J2(?z )u ] u(° 2 ° 2 \U ] Ue V Vq D 12 Ji<Jz ) jj^U 

V °2 cr z J 



1,1 



J\°x Vqd12 
VqdI-2 J2&z 



2 2 

O2 cr z 



J Jl°x Vqd12 
V ^QDl2 Jl&z 



[si(-t)] 2) 2. 



J 2,2 



Taking into account that [X2(t)]i 1, p2(£)]2 2 are even functions of £, Eq. (23) follows. 

I — 



(A2) 



APPENDIX B: PROPAGATION ALGORITHM 

Let H{t) = E a H a (t) + Za(H° aC + H Ca ) + Hc(t) be 
the one-particle Hamiltonian of a system which consists 
of a = 1,2,. ..,7V electrodes in contact with a central 
region C. We assume that the time dependence of 



H a (t) 



K1 

H° n 



U a ,i(t) 
U ad (t) 



Hi 



U a (t) 



is a uniform spin-dependent and time-dependent shift 
while the time-dependence of He has no restrictions. 
We denote with \ip a ) the projection of a generic wave- 
function \ip) on electrode a and with \tpc) the projection 
of 1^) onto region C. The time-dependent Schrodinger 
equation reads 

i^hM*)> = H a {t)\i, a {t)) + H° aC \^ c (t)), (B2) 

i^c{t)) = H c (tMc{t)) + J2 H cMt))- (B3) 



Performing the gauge transformation 

\1> a (t))=exp[-i f&rU a {r)] \(/> a {t)), (B4) 
Jo 



and \ip c (t)) = |0cW), Eqs. (|B2][B3| become 
. d 



i-\K{t)) = H°Mt)) + H aC mc(t)), (B5) 
(Bl) ^Mc{t)) = ifcWI^c^ + E^WI^W)'^ 6 ) 



with H c ait) = H° Ca ex.p[-i J^drU a (r)] and H aC (t) = 
[Hca(i)y. The effect of the gauge transformation is to 
transfer the time dependence from the Hamiltonian de- 
scribing the bulk electrodes to the Hamiltonian describ- 
ing the contacts between the electrodes and region C. 
The gauge-transformed Schrodinger equation is used to 
calculate the time evolved state \4>{t m = mAt)) = \<j)^ m ^) 
by using the Cay ley method 

(1 + i5H g m) )\^ m ^) = (1 - iSH g m) )\^), (B7) 
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where 5 = At/2, = \[H g {t m+1 ) + H g (t m )}, and 

H g {t) = E a H a +J2 a (H aC it)+H Ca (t))+Hc(t) is the 
gauge-transformed Hamiltonian. The interface Hamilto- 
nian H a c is spin-diagonal provided region C includes 
the first few atomic layers of electrode a. In this case the 



projection of Eq. (B7) onto different subregions leads to 



a close recursive relation for the amplitudes |</>[^) of the 
wave-function in region C (the steps are similar to those 
of Ref. 



1 'A£7~( m ) 

| ^ ro+1)) = i-ZdM ^ )} + |5(m)) _ |M(m)) (B8) 

i + idHy 

where the source term \SW) and the memory term 
|M( m )) read 



In the above equations we have used the following defi- 
nitions 



„(m) 



exp[-z J £r7l+1 dr U a [r)\ + exp[-i f* m dr U a (r)} 

2 ' 
(Bll) 



n(m) _ ttO (1 - iSH^)™ o 



H 



(m) 
eff 



H 



(m) 



Q 



(0)^(m)_ 



(B12) 



(B13) 



|S (m) ) 



1 + WJf 



2 ^ V^^(m)rrO (1 ~ ^(5iJ^) m (Q) 



(m) 

eff a 



(B9) 



|M (m) ) 



m—1 



i-r-zo.rz eff q/ j=0 



EE*W+^ +1) ) 



x 2; 



(-^(I^Vl^" 1 "^)- (BIO) 



The recursive relation in Eq. (B8) is written in terms 



of matrices and vectors with the same dimension of the 
central region, i.e., the infinitely large electrodes have 
been embedded in an effective equation of finite dimen- 
sion. We defer the reader to Ref. [25] for the description of 
how to calculate the matrices and the source term 

| S (m)^ 
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